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A  METHOD  OF  CALCULATING  SUBLIMATION  RATE  FOR 
THERMIONIC  CATHODES 

Guo  Wenxiang,  Institute  of  Electronics,  Chinese  Academy  of 
Sciences 

[Abstract]  An  improved  method  of  measuring  cathode 
sublimation  rates  using  a  crystal  oscillator  is  proposed;  in  the 
method,  the  distance  between  the  crystal  (substrate)  and  cathode 
(source)  is  smaller  than  is  usual.  Thus,  more  sublimation 
products  from  the  cathode  (source)  can  be  received  by  the  crystal 
and  consequently  the  method  can  serve  in  measuring  lower 
sublimation  rates.  Formulas  for  calculating  the  sublimation  rate 
are  derived  and  their  precision  is  stated.  In  the  paper,  a 
method  of  preventing  the  resonant  frequency  change  caused  by 
temperature  changes  is  also  given.  The  experimental  results  show 
that  1Hz  of  shift  in  the  resonant  frequency  can  be  detected  when 
the  oscillation  frequency  is  5MHz. 

I .  Introduction 

By  using  the  crystal  oscillator  method  in  measuring 
sublimation  rates  of  cathodes,  generally  the  cathode  is  regarded 
a  point  sublimation  source  [1,2]?  this  is  possibly  because  of 
facilitating  the  calculation  of  the  sublimation  rate.  In  other 
reports,  such  as  in  reference  [3],  the  cathode  is  considered  as  a 
flat-element  sublimation  source;  the  sublimate  receiver  is 
considered  as  a  flat  element  at  a  location  along  the  normal 
direction  of  the  cathode  surface. 
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The  article  adopts  an  effective  measure  for  overcoming  the 
effect  of  measuring  the  sublimation  rate  by  the  oscillation 
frequency  drift  (crystal  resonant  frequency  drift)  of  the  crystal 
because  of  temperature  change.  This  not  only  enhances  the 
precision  in  measuring  the  oscillation  frequency  change  (crystal 
resonant  frequency  change)  due  to  the  sublimate;  in  addition,  it 
allows  a  closer  distance  between  the  crystal  and  cathode.  That 
is,  in  a  situation  of  not  changing  the  area  of  the  crystal 
receiving  the  sublimate,  the  quantity  of  received  sublimate  per 
unit  time  is  increased.  In  this  article  the  cathode  is 
considered  as  a  flat-element  sublimation  source.  Under  the 
conditions  of  greater  area  ratio  in  receiving  the  sublimate, 
equations  and  inaccuracies  in  calculating  the  cathode  sublimation 
rate  are  derived.  In  addition,  the  actual  application  situations 
are  introduced. 

II.  Premise  Conditions  of  Calculation 

It  is  assumed  that  an  equilibrium  state  is  received  between 
the  molecules  just  sublimated  from  the  cathode,  and  the 
sublimated  molecules  in  the  outside;  at  that  instant,  the 
sublimation  rate  is  equal  to  the  condensation  rate.  If  the 
sublimated  molecules  are  not  returned,  then  the  sublimation  rate 
(condensation  rate)  can  be  measured  [4].  Numerically,  this  is 
the  collision  rate  of  the  sublimated  molecules  against  the 
cathode  in  the  equilibrium  state.  Computing  the  collision  rate 
of  the  flat  element  by  gas  molecules  in  the  equilibrium  state  is 
the  general  method  of  calculating  the  sublimation  rate  of  flat- 
element  sublimation  sources. 

In  the  following,  particular  situations  of  cathode 
sublimation  are  analyzed.  Considering  that  the  upper  limit  (the 
maximum  rate)  of  the  cathode  sublimation  rate  is  smaller  than 
1016atoms/cm2/s  (if  barium  atoms  are  sublimated,  this  upper  limit 
of  cathode  rate  corresponds  to  a  value  smaller  than 
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2.3’10"6g/cm2/s) ,  then  the  corresponding  vapor  pressure  of  the 
cathode  sublimate  is  lower  than  10“4torr  [5].  If  the  effective 
diameter  of  the  sublimated  molecules  is  5*10“8cm,  and  the 
temperature  is  higher  than  1000K,  then  the  mean  free  path  length 
is  greater  than  90cm.  If  the  cathode  size  is  the  conventional 
cathode  to  be  measured,  and  the  cathode  is  in  a  vacuum,  then  if 
the  cathode  sublimation  is  considered  as  molecules  in  the 
equilibrium  escaping  into  the  vacuum  [6]  through  small  holes  from 
the  vessel,  then  because  of  the  condition  that  the  mean  free  path 
length  is  considerably  greater  than  the  cathode  dimension,  the 
equilibrium  state  of  molecules  in  the  vessel  can  be  considered 
such  that  there  is  no  effect  of  molecules  escaping 
(notwithstanding  that  the  total  number  of  molecules  in  the  vessel 
is  slowly  decreasing)  and  the  molecules  escaping  in  the  form  of 
effusion.  Therefore,  the  calculation  of  the  sublimation  rate  can 
be  considered  as  the  calculation  of  collision  number  per  unit 
area  of  the  cathode  surface  by  molecules  per  unit  time  at  the 
corresponding  vapor  pressure. 

III.  Calculation  of  Cathode  Sublimation  Rate 

Assume  that  the  area  receiving  the  sublimate  is  a  circle 
with  diameter  2a,  and  that  the  sublimation  area  of  the  cathode  is 
the  area  element  ds,  and  two  circles  are  perpendicular  to  the 
connecting  line  oo'  of  the  circles'  center  (Fig.  1).  Then  the 
collision  number  (refer  to  Fig.  1)  to  ds  per  unit  time  for 
molecules  in  the  three-dimensional  angle  (the  three-dimensional 
angle  of  a  cone  with  vertex  angle  of  20q)  of  the  corresponding 
circle  area  receiving  the  sublimate. 

From  the  theory  of  molecular  motion,  the  number  of  molecules 
colliding  at  ds  per  unit  time  for  a  three-dimensional  angle 
element  du  with  angle  9  formed  with  oo'  is 

mi  cot  6  it,  (  1  ) 

In  the  equation,  Z  is  the  average  speed  of  molecules  and  n  is  the 
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molecule  number  per  unit  volume.  By  establishing  a  spherical 
surface  coordinates  with  point  o  as  the  origin,  then  du>>sin0d6dt. 
Then  the  number  of  molecules  per  unit  time  colliding  at  ds  of  the 
origin  o  of  the  three-dimensional  angle  of  the  corresponding  cone 
with  vertex  angle  26q  is? 


tin  5  co*  848 


—  »n,8t4t. 


(2) 


In  the  equation,  nv/<  is  the  overall  molecules  in  the  direction 

from  the  upper  side  of  ds,  colliding  at  a  unit  area  ds  per  unit 

p 

time.  This  can  also  be  expressed  as  ^/~Tnm{T  (P  is  the 

vapor  pressure,  k  is  the  Boltzmann  constant,  m  is  molecular 
weight,  and  T  is  absolute  temperature).  Numerically,  np/4  is 
equal  to  the  sublimation  rate  at  ds.  Then  if  N0  is  known,  the 
sublimation  rate  can  also  be  known;  however,  N0  can  be  measured 
experimentally . 


Fig.  1  Fig.  2.  Schematic  diagram  of 

experimental  tube  for  measuring 

cathode  sublimation  rate 

KEY:  1  -  Crystal  2  -  Screen  plate 

3  -  Movable  back  plate 

4  -  Cathode 


The  schematic  diagram  for  using  the  crystal  oscillator  for 
measuring  the  cathode  sublimation  rate  is  Fig.  2.  In  the  figure, 
the  purpose  of  the  screen  plate  is  to  control  the  area  in  which 
the  sublimate  is  received  by  the  crystal.  The  purpose  of  the 
movable  back  plate  is  to  open  or  cut  off  the  passage  of  the 
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sublimate  at  any  time  in  order  to  precisely  control  the 
sublimation  time.  D  is  the  diameter  of  the  area  receiving 
sublimate  by  the  crystal,  d  is  the  cathode  diameter,  a  is  the 
radius  of  the  screen  plate  hole,  and  b  is  the  distance  between 
the  screen  plate  and  the  cathode. 


By  utilizing  the  triangular  relationship  shown  in  Fig.  2, 
the  molecule  number  N0  is  replaced  by  the  corresponding  mass  Mq/ 
the  sublimation  rate  is  expressed  in  v,  and  only  ds  at  the  center 
of  the  cathode  surface  is  considered,  then  from  Eg.  (2)  we 
obtain : 


r_  */.(«»  +  *’) 

m'J, 


(3) 


Considering  that  the  actual  sublimation  time  of  measurement 
is  t,  the  cathode  is  a  cylindrical  terminal  surface  emission 
cathode  with  area  s  (diameter  d),  then  the  average  sublimation 
rate  of  the  time  period  t  is  as  in  the  following  approximate 
relation: 


r  _  «.•(.’  +  *’) 
- - , 


(4) 


In  the  equation,  M*  is  the  mass  of  the  sublimate  received  by  the 
crystal  from  actual  sublimation  of  a  limited-size  cathode. 


After  the  crystal  received  the  sublimate,  the  variation 
(decreasing  value)  df  of  the  oscillation  frequency  can  be 
expressed  as  [ 6 , 7 ] :  . 

(5) 

In  the  equation,  F  is  the  area  in  which  sublimate  was  received  by 
the  crystal;  dm  is  the  mass  of  the  sublimate  received  by  the 
crystal;  Cf  is  the  quality  determination  sensitivity,  which  can 
be  experimentally  determined  according  to  the  actual  application 
situation.  Under  the  above-mentioned  conditions  of  measuring  the 
cathode  sublimation  rate,  this  relation  is  derived.  Since  the 
adherent  thickness  of  the  sublimate  is  within  the  range  of 
hundreds  of  angstroms  or  even  ten  thousand  angstroms,  df  and  dm 
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have  a  better  linear  relationship  [1,7,8].  Therefore,  from 
Eq.  (4)  and  utilizing  Eq.  (5),  we  obtain: 

,  -  <«  ± g)/PVA/ 

'  .V,  \7/  "T*  (6) 

in  the  equation,  Af  indicates  the  decreasing  value  of  oscillation 
frequency  after  the  crystal  received  the  sublimate. 

IV.  Inaccuracies  of  vt 

After  the  actual  cathode  is  used  to  replace  a  flat  element, 
the  accuracy  thus  created  can  be  estimated  as  follows. 

As  shown  in  Fig.  3,  if  the  area  in  which  the  sublimate  is 
received  is  s0  (2a  is  the  diameter),  the  cathode  area  is  s  (d  is 
the  diameter),  and  two  circles  are  perpendicular  to  the 
connecting  line  oo'  of  the  circles'  center.  Passing  through  the 
cathode  edge  at  a  point  o1#  o.o.'l oo,  ,  by  making  go',  —  o\h,  , 
then  the  area  element  of  the  point  o^  is  for  the  sublimation 
quantity  .  (refer  to  Eq.  (2))  of  area  S!  (the 

diameter  is  2 (a  -  r),  where  r  is  the  cathode  radius),  for  the 
area  element  of  the  point  Oj^.  Actually,  the  sublimation  quantity 
No?|Ni  of  an  area  element  of  the  point  o^  that  sq  receives, 
because  other  than  s^,  the  area  (sq  -  s^)  already  receives  the 
sublimate  of  the  area  element  of  point  oj_. 

Assume  that  the  sublimation  quantity  N0^Nt  —  ~ua,6^tt 
refer  to  Eq.  (2)),  for  sq  by  area  element  at  point  o,  then  Nq>N'. 
The  following  can  be  proved:  refer  to  Fig.  3,  make 

/0;  -  ,,  . ,,  .  Obviously,  the  sublimation  quantity 

N2=Hq  to  s2  by  area  element  at  point  oj_.  Since  whether  Sq  or  s2, 
for  the  sublimation  of  area  element  at  point  o^,  the  sublimation 
quantities  are  the  same  for  the  area  (the  portion  marked  by 
hatching  in  the  figure)  between  fh,  however,  the  included  angle 
(e>e2)  between  the  partial  cone  and  OiOi  corresponding  to  area  I 
between  ef;  the  included  angle  g<e2  between  the  partial  cone  and 
o^o(  corresponding  to  area  II  between  hi.  In  addition,  since  the 
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three-dimensional  angle  of  the  area  I  to  o^  is  smaller  than  the 
three-dimensional  angle  of  II  to  o^ ,  then  according  to  Eq.  (1) 
NqcNq  can  be  obtained.  Therefore,  N1<Nq<N0. 


Fig.  3  Fig.  4 

It  can  be  proved,  with  the  same  line  of  reasoning  as  above, 
that  to  any  point  on  between  oo^ ,  there  obtains  the  relation 
Nn<N0<N0  (»o  is  the  sublimation  quantity  of  area  s0  for  the  area 
element  of  point  on;  Nn  is  the  sublimation  quantity  of  area  sn  of 
area  element  at  point  on.  The  diameter  of  sn  is  2 (a  -  rn). 
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rn=00n)-  Since  N1<Nn,  the  boundary  of  inaccuracy  is  (Ng-N^ )/Ng 
for  N0  of  the  sublimation  of  area  elements  of  various  points  to 

Sg  . 


If  we  apply  Eq.  (4)  or  Eq.  (6)  to  calculate  the  sublimation 
rate,  actually  some  number  greater  than  but  smaller  than  Ng  is 
considered  as  Ng ,  therefore  the  boundary  of  inaccuracy  of  v^-  can 
be  expressed  as: 

Mo  -  M,  _  No  -  N,  _  j _ 1 _ 

M,  N,  a 3  , _ o'M  (  7  ) 

«J+*J  („’  +  **)(.-  r)* 

In  the  following,  consideration  can  be  applied  to  shrink  the 
boundary  of  inaccuracy.  Refer  to  Fig.  4,  by  using  point  o^  as 
the  center  of  the  circle,  and  o^l  as  the  radius  (1  is  the 
intersection  between  o^k^  and  the  circle  by  using  o'  as  the 
center  of  the  circle  and  a  as  the  radius,  among  o^k,1//o/k# ) . 

Then  based  on  Eq.  (2)  or  Eq.  (3),  the  sublimation  quantity  of 
area  jlno^pinj  for  area  element  of  point  O!  (refer  to  Fig.  3) 
should  be 


1  tiv  a1  —  r J  1 

2  4  (a1  —  r1)  +  b1  J  ‘  ;  similarly,  the  sublimation  quantity  of  area 

o^nhpo^  for  the  area  element  of  point  is  JL  [  (a  —  ry  i 

2  4  17.  -  ry  +  y  J  dt 

Thus,  it  is  found  that  the  sublimation  quantity  N^  of  area 
jlnhpmj  for  the  area  element  of  point  o^  is 


-L  ™  [  **  r>  +  (a  -  rV  ]  . 

2  4  I  (V  _  (a_r)j  +  ^i  J  .  Therefore,  with  the  above- 

mentioned  analysis,  it  is  found  that  the  boundary  of  inaccuracy 
6vt  is  further  shrunk: 


u, -  - 1  - 

*  N,  2a1 


(.  -  r)> 


Consideration  is  given  to  the  fact  that,  during  the  actual 
measurement,  0<r<a,  therefore,  a  and  b  are  determined,  6v^ 
increases  simply-harmonically  with  an  increase  in  r.  As  an 
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example,  if  we  take  a=0.40cm  and  b=0.50cm,  when  r  is, 
respectively,  0.15cm,  0.10cm,  and  0.05cm,  according  to  Eq.  (a), 
ovt  is,  respectively,  29%,  18%,  and  8%;  if  we  take  a=b=0.50cm,  and 
r  is  still,  respectively,  0.15cm,  0.10cm,  and  0.05cm,  then 
according  to  Eq.  (8)  6v^  is  obtained,  respectively,  equal  to  19%, 
12%,  and  5%.  (The  boundary  of  the  actual  sublimation  rate 
greater  than  v^  does  not  exceed  <5v^ . ) 

V.  Applications 

The  above  analysis  adaptable  to  the  situation  in  which  the 
crystal  is  at  a  relatively  close  distance  from  the  cathode;  at 
the  same  time,  the  effect  on  crystal  temperature  of  thermal 
radiation  from  the  cathode  is  also  intensified.  The  rise  in 
crystal  temperature  will  cause  drifting  of  the  oscillation 
frequency:  this  is  a  very  key  problem  in  measuring  the  effect  of 
trace  amounts  of  sublimate.  According  to  reports  in  reference 
[2],  within  the  range  of  -40  to  -t-90°C  by  cutting  crystals  for  AT 
with  higher  thermal  stability,  the  frequency  drift  is 
approximately  plus  or  minus  40Hz.  When  the  temperature  exceeds 
100° C,  the  frequency  drift  rises  rapidly  with  increase  in 
temperature.  According  to  the  ordinary  sublimation  rate  of  the 
cathode  and  the  geometry  of  the  assumed  experimental  tube,  an 
estimate  is  made  that  the  sublimation  time  required  for  the 
frequency  to  be  lowered  to  40Hz  is  approximately  several  hours. 

In  other  words,  if  we  desire  that  the  sublimation  time  not  to  be 
in  excess  of  several  hours,  the  frequency  change  smaller  than 
40Hz  should  be  detected.  To  avoid  the  frequency  drift  affecting 
measurements  due  to  temperature  change,  in  some  cases  the  crystal 
is  to  be  placed  further  from  the  cathode  (approximately  6  to 
14inches  as  the  distance)  and  apply  cooling  to  the  crystal 
(maintain  the  crystal  temperature  lower  than  100°C)  [2].  In  some 
cases,  two  identical  crystals  are  used  to  operate  in  identical 
environments.  In  one  crystal,  the  sublimate  is  received;  in 
another  crystal,  by  using  a  screen  as  a  barrier  to  most  heat 
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radiation  past  through,  not  to  receive  the  cathode  sublimate.  By 
measuring  the  frequency  difference  of  the  oscillation  frequencies 
of  these  two  crystals,  this  is  the  frequency  decrease  value  [1] 
caused  by  the  sublimate  due  to  a  frequency  difference  of  these 
two  crystals.  In  the  article,  frequency  measurements  are 
conducted  before  and  after  the  sublimation  to  keep  the  isothermal 
crystal  at  the  same  temperature  (at  a  certain  temperature  nearing 
the  point  of  the  zero  frequency  temperature  coefficient),  the 
frequency  difference  is  taken  as  the  value  of  frequency  reduction 
caused  by  the  sublimate  thus  received.  Fig.  5  shows  the 
frequency  temperature  curve  actually  measured  at  different 
temperatures  for  an  isothermal  crystal.  The  experiment  was 
conducted  on  a  chip  of  crystal  sealed  in  an  experimental  tube 
from  which  air  had  been  expelled.  The  marks  in  the  figure  show 
full  data  measured  on  the  same  day  (restart  of  the  crystal 
oscillator).  From  Fig.  5,  if  the  isothermal  temperature  during 
measurement  is  within  the  range  (approximately  between  64.5  and 
68.5 °C  in  the  fig-  re)  of  plus  or  minus  2°C  for  the  measured 
isothermal  temperature  at  the  zero  frequency  temperature 
coefficient  point,  the  fluctuation  of  the  crystal  isothermal 
temperature  is  controlled  not  to  exceed  plus  or  minus  0.5°C,  then 
the  range  of  frequency  variation  can  be  smaller  than  0.5Hz  caused 
by  the  fluctuation  of  isothermal  temperature  and  frequency 
reappearance  (such  as  temperature  increase  and  decrease  for  the 
crystal  and  restarting  of  the  crystal  oscillator).  Consideration 
is  given  to  the  results  of  actual  measurement  on  frequency 
stability  of  the  crystal  oscillator  when  the  oscillation 
frequency  is  5MHz;  the  range  of  frequency  variation  can  be 
smaller  than  0.2Hz/24h  (not  considering  the  stability  of  the 
frequency  measurement  instruments),  after  another  consideration 
is  given  that  the  crystal  received  the  cathode  sublimate  on  the 
frequency-time  relationship  at  a  certain  temperature  near  the 
zero  frequency  temperature  coefficient  point,  the  experiment 
measured  that  the  frequency  approximately  was  reduced  by  0.5Hz 
after  24h?  however,  when  several  hours  were  actually  required  in 
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measuring  the  sublimation  rate  at  a  certain  temperature. 
Therefore,  the  detection  change  of  frequency  precise  to  1Hz  (when 
the  oscillation  frequency  is  5Mhz)  is  possible.  Now,  we  consider 
cf  [ 6,7],  in  Eq.  (5)  we  know  that  the  variation  of  frequency  for 
1Hz  corresponds  to  barium  with  a  single  atomic  layer. 


(Hr) 

5000671.0  ■ 

5000670  0  ■  « 

5000669.0  - 
50006680  - 
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—i - 1 - 1— 
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Fig.  5.  Crystal  oscillator  oscillation  frequency 
versus  crystal  temperature 


VI.  Conclusions 


1.  In  this  method,  the  sublimate  is  received  within  a 
relatively  short  distance  for  a  relatively  large  area,  therefore 
the  efficiency  of  receiving  the  sublimate  is  increased,  capable 
of  effectively  measuring  relatively  low  sublimation  rates.  By 
utilizing  Eqs.  (6)  or  (4),  the  sublimation  rate  can  be 
conveniently  calculated  for  a  limited  cathode  area.  The  boundary 
of  imprecision  for  the  sublimation  rate  calculation  formula  is 
determined  by  Eq.  (8)  (in  the  derivation,  it  was  assumed  that  the 
sublimation  rate  is  the  same  at  various  portions  of  the  cathode 
surface).  During  applications,  requirements  can  be  based  to 
select  the  values  a,  b,  and  r. 

2.  In  this  article,  measures  are  taken  to  overcome  the 
crystal  oscillation  drift  caused  by  thermal  radiation  of  the 
cathode,  thus  basically  eliminating  the  effect  on  measurements 
due  to  frequency  drift  caused  by  the  crystal  temperature  rise  or 
the  crystal  temperature  change,  thus  enhancing  the  precision  of 
frequency  change  due  to  measurement  of  sublimate.  As  indicated 
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by  the  experiments,  when  the  oscillation  frequency  is  5MHz,  the 
measured  frequency  change  can  be  precise  to  1Hz. 

The  article  was  first  received  on  August  11,  1981. 

The  revised  article  was  received  on  June  27,  1983. 
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